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Abstract. We have recently started to investigate 2D arrays of confocal lens 
pairs. Miniaturization of the lens pairs can make the array behave ray-optically 
like a homogeneous medium. Here we generalize the geometry of the lens pairs. 
These generalisations include a sideways shift of the lens centres and a change in 
the orientation of both lenses in each pair. We investigate the basic ray optics of 
the resulting arrays, and illustrate these with movies rendered using ray-tracing 
software. We suggest that confocal lenslet arrays could be used to realize ray- 
optically some recent metamaterials concepts such as the coordinate-transform 
design paradigm. 
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1. Introduction 

The invention and experimental realization of metamaterials with highly unusual 
optical properties [H El |3] such as negative refraction and resonant amplification of 
evanescent waves has opened up one of the most exciting current research areas in 
optics. Metamaterials have led to many new concepts, for example superlenses that 
image not only travelling waves but also evanescent waves [4] , and a new optical-design 
paradigm based on coordinate transforms [5l |6] . 

The manufacture of metamaterials is still exceedingly difficult [7 . Even the 
simplest metamaterials for near- visible wavelengths (e.g. [8 ) are currently not very 
big (interference lithography can structure areas of order Icm^ |7]) and difficult to 
manufacture. Metamaterials are only now beginning to work over a broad band 
of frequencies [9]. Our work is motivated by the idea of using standard optical 
components - lenslet (or microlens) arrays - to enable many (but not all) of the 
unusual ray-optical properties of metamaterials to be "experienced" by approximating 
these metamaterial properties ray-optically on an everyday scale (~10cm) and across 
the entire visible spectrum. Perhaps most excitingly, our ray-optical materials enable 
light-ray direction changes without wave-optical analog [TOl. Their realization could 
therefore open up entirely new concepts. 

Like with metamaterials, our starting point is negative refraction, specifically the 
recent realization [11^ that a sheet formed by two identical, parallel, lenslet arrays 
(LAs) with a common focal plane (so the lenslet arrays are confocal) [12] |13] acts 



Generalized refraction using lenslet arrays 



2 



ray-optically like the interface between optical media with opposite refractive indices. 
(Note that wave-optically such a sheet does not act like the interface between optical 
media with opposite refractive indices. There is, for example, no negative group 
velocity [14] and no amplification of evanescent waves required for the sub- wavelength 
imaging properties of super lenses [4] and hyper lenses [15l[8].) A sheet of two confocal 
LAs (CLAs) with different focal lengths approximates the interface between optical 
media with a refractive-index ratio given by minus the focal- length ratio [TJ. The 
refractive-index ratio can be positive or negative. This holds only for the part of the 
light that passes through corresponding lenslets [16] . CLAs can be designed such that 
other light is absorbed; if it is not absorbed, it leads to "ghost images" [16 . 

Here we generalize CLAs. Simple modifications (such as a sideways displacement 
of the centres of corresponding lenslets in the two LAs) offer additional design 
parameters that will enable unprecedented control over light rays, while retaining 
the important property of the structure acting macroscopically like a homogeneous 
material. Generalized CLAs can then do significantly more than simply act like the 
interface between different refractive indices, and this paper provides a concise atlas 
of distorted views through different individual CLAs. 

We speculate that generalized CLAs can offer control over light-ray propagation 
similar to that afforded by metamaterials. This could enable, for example, a ray- 
optical approximation of meta-material concepts such as the coordinate-transform 
design paradigm [6 . Because of this and other analogies with metamaterials we call 
these generalized CLAs meta materials for rays (METATOYs) [TO] , 

2. Generalization of confocal lenslet arrays 

A CLA sheet consists of two LAs, one with focal length /i, the other with focal length 
/2, arranged so that they share a common focal plane. 

CLAs with /i = /2 have long been used. The first application was 3D imaging, 
initially for "integral photography" [17 , more recently for 3D displays (see Ref. [T8] 
and references therein). Related to imaging is the moire magnifier [19 , which is 
based on "mis- aligned" CLAs. Without making any connection with metamaterials, 
what can be argued to be the analogs of a superlens [H [20] and a hyperlens 
[m [i [E] have already been built using CLAs [H [13 [H |23] . However, these do 
not offer the amplification of evanescent waves and corresponding sub-wavelength- 
resolution-imaging properties afforded by their metamaterial counterparts, and they 
are therefore arguably much less interesting. Nevertheless, all of these are examples 
of the possibilities afforded by CLAs. 

An important property of the CLAs we consider here is that the change of ray 
direction is independent of the precise point where a light ray hits the array. If the 
offset of the light rays on transmission through the CLAs is small (as is the case if the 
individual lenslets are small [11 ), CLAs therefore have a homogeneous appearance. 
This property can also be realized with arrays of prisms ^24| |25| . In CLAs it is due to 
the following simple argument. Two parallel light rays, after hitting the first lenslet, 
are focussed so that they pass through the same point of the lenslet's back focal plane. 
But as this plane is also the second lenslet's front focal plane, the two light rays are 
parallel again after passing through the second lenslet. This is not normally the case, 
like for example in a single LA or non-confocal pairs of LAs. 

To be able to realize metamaterials concepts that generalize refraction, CLAs 
need to be generalized. The following generalisations keep corresponding lenslets in 
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Figure 1. Geometry and ray optics of generalized confocal lenslet arrays (CLAs). 
(a) Schematic of two-dimensional generalized CLAs. An individual lens (focal 
length /i) in the left lenslet array focusses all light rays that impinge on it with 
the same angle of incidence a.i in a point in the focal plane F. These rays leave 
the corresponding lens (focal length /2) in the right array with an exit angle ol2- 
The light rays that pass through the centre of each lens are shown as a dotted 
line. Both lenses, and with them their common focal plane, F, are inclined with 
respect to the plane of the CLAs by an angle 6. In addition, the optical axes of 
the lenses are displaced by a distance d. (b) Ray-tracing image of a small part 
of three-dimensional generalized CLAs. The picture shows the top edge of the 
two lenslet arrays, rendered as transparent with a blue tinge, in front of a light- 
brown/dark-brown checker-board floor. Each lenslet array consists of identical 
lenslets, but the lenslets that make up the two arrays are different. This image 
was rendered using the ray-tracing software POV-Ray [26) . 

the two arrays confocal and therefore the appearance of the CLAs homogeneous, 
while introducing additional parameters that affect the way the light-ray direction is 
changed: 

• the centres of corresponding lenslets can be offset with respect to each other; 

• each pair of corresponding lenslets, and with them their common focal plane, can 
be inclined; 

• the lenslets can be made elliptical, whereby the focal lengths in the directions 
of the axes of the ellipses of corresponding lenslets have to add up to the same 
value, namely the separation between the lenses. 

Fig. [T] shows examples of such generalized CLAs. Manufacturing such generalized 
CLAs will require a significant extension of current techniques for manufacturing 
complex microlens arrays [27^ 28 . 

From now on, we will usually drop the adjective "generalized" and simply refer 
to CLAs. In the following sections we discuss CLAs in two and three dimensions. 

3. Two-dimensional case: cylindrical CLAs 

It is instructive to start by considering the two-dimensional (2D) case. This case also 
describes transmission through cylindrical CLAs (more precisely, confocal arrays of 
cylindrical lenslets), specifically the projection into a plane normal to the cylindrical- 
lenslet axes. 

Figure [ija) shows a cross-section through CLAs in two dimensions. The two 
central lenslets have focal lengths /i and /2, and they are separated by the sum of 
their focal lengths, /i + /2. The two lenslets' optical axes are parallel; that of the 
second lenslet is shifted with respect to that of the first by a distance d. In addition. 
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both lenslets, and with them their common focal plane, F, are inclined with respect 
to the vertical (the plane of the CLA sheet) by an angle 0. This means the optical 
axes are inclined by the same angle with respect to the sheet normal. 

Light rays hitting the left lenslet at an angle ai with respect to the sheet normal 
(the horizontal in Fig.jlja)) have an angle (ai —6) with respect to the lenslet's optical 
axis. They therefore pass through the common focal plane a distance 

c = fi tan(ai — 0) (1) 

from the left lenslet's optical axis. Provided they then pass through the corresponding 
lenslet on the right, such rays leave at an angle (0^2 — 0) with respect to the right 
lenslet's optical axis such that 

tan(a2-6>) = (d-c)//2. (2) 

a2 is the angle of refraction: the angle with respect to the sheet normal at which the 
rays leave the sheet. Eliminating c from these equations gives 

/i tan(ai — 0) = d — f2 tan(a2 — 0). (3) 

With the dimensionless quantities 

. d 



and 



, , (4) 



Jl 



this becomes 

tan(ai — 0) = S -\- 7^tan(a2 — 0). (6) 

Equation (|6| is one of the key results of this paper. It describes a generalized "law of 
refraction" for CLAs in two dimensions, or cylindrical CLAs. In the following section 
we generalize this to the three-dimensional case. 



4. Three-dimensional case: elliptical CLAs 

We start our generalisation of the 3D case with two confocal spherical lenses with a 
common optical axis, the z axis. Such a pair of confocal lenses has the important 
property that the direction of the refracted ray depends only on the direction of the 
incident ray and is independent of the exact position where the ray hits the first 
lens [11 . In this section we discuss generalisations of the lens pair that retain this 
property, and we consider the two lenses to be corresponding lenslets in CLAs. Fig. 
[TJb) shows a ray-tracing image of such generalized CLAs. It was created using the 
free, open-source, ray-tracing program POV-Ray [26 , which we also used throughout 
the remainder of this section to illustrate the view through examples of generalized 
CLAs. 

As the first step of our generalisation we replace the spherical lenses with elliptical 
(or astigmatic) lenses. We consider the case of elliptical lenses whose major and minor 
axes respectively point in the x and y directions or vice versa. (The restriction of the 
major and minor axes to the x and y directions will be generalized later.) The effect 
of an elliptical lens on a light beam can be understood in the following, idealized. 
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Figure 2. View through confocal lenslet arrays with different values of rjx- Each 
frame shows the simulated view of a chess piece seen through CLAs. The CLAs 
are positioned in the z = plane, the chess piece is positioned at = +10, the 
camera is positioned at Zc = —20. In this figure (and all subsequent figures), all 
lengths are given in units of the fioor-tile length. The CLAs consist of 200 x 200 
confocal lenslet pairs, covering a square area of side length 1. The separation 
between corresponding lenslets is 0.05; all lenslets are convex-convex with a centre 
thickness of 0.0005. (We restrict ourselves throughout this paper to simulating 
convex-convex elliptical lenslets, which we describe as the intersection between 
two similar, but displaced, ellipsoids.) All views were simulated with POV-Ray 
[26| . In subsequent figures we use the case rjx = —1/4, r]y = —1 as the starting 
position for further generalisation. 



wa}||] An elliptical lens is approximately equivalent to two orthogonal cylindrical 
lenses aligned with the elliptical lens's major and minor axes. We consider the effect of 
the elliptical lens on a transmitted ray bundle in terms of the orthographic projections 
into the x-z and y-z planes. In each projection we only consider the effect of the 
cylindrical lens whose corresponding cylinder axis is normal to the projection plane, 
that is, the lens which has curved surfaces in the projection; we ignore the effect of 
the other cylindrical lens. The lens then has two different focal lengths in the two 
projections, fx in the x-z projection and fy in the y-z projection [29 . 

To retain the position independence of the ray-direction change, the two lenses 
have to be confocal in both projections. In other words, the focal lengths in the x-z 
projection of the first and second lens have to add up to the lens separation, and the 
same has to be true in the y-z projection, so 

fxl + fx2 = fyl + fy2' (7) 

In analogy to the previous definition, we describe the focal-length ratios in the x-z 
and y-z projections by the dimensionless quantities 

fx2 fy2 /Q\ 

r]x = -^^ Vy = -j-' (8) 

Jxl Jyl 

i Note that the standard geometrical treatment of the effect of spherical lenses on light rays, which 
we have used in Ref. JT, and in the previous section to derive Eqn JgI, is also idealized. 



Generalized refraction using lenslet arrays 



6 




1 






































z,=28 



z,=38 



Figure 3. Dependence of the apparent stretching of an object seen through 
CLAs on the distance between the CLAs and the object. Here the object is a chess 
piece placed at various distances zj^ behind the CLAs. The frames shown here are 
taken from a movie showing the view through the CLAs as the chess piece moves 
between = —2 to Zk =38. The movie is contained in the supporting onhne 
material (MPEG-4, 184kB). This simulation was performed with rjx = —1/4 and 
all other parameters like in Fig. [2] 



Fig. |2] shows a chess piece as seen through different CLAs, ah with r]y = —I but with 
different values of r]x- Most of the visible part of the chess piece (everything apart 
from the plinth) is positioned in a plane parallel to the CLAs and appears squashed 
or stretched in the horizontal {x) direction. 

Fig. [2] and all subsequent figures are calculated for CLAs comprising pairs of 
confocal lenslet pairs with a relatively small (1:100) ratio between lenslet diameter 
and distance between the lenslets. This ratio allows all lenslets to be relatively thin 
and used close to their optical axis, which leads to a respectable image quality. On the 
other hand, it severely limits the field of view of the CLAs (approximately to a cone 
with a cone angle of order tan~^ (1/100)). Consequently, all figures are calculated as if 
the scene was photographed with a lens with a long focal length, and correspondingly 
small field of view. More careful lens design should enable a significant increase in the 
CLA's field of view without compromising optical quality. 

Fig. [3] demonstrates that the squashing and stretching shown in Fig. [2] depends on 
the distance between the CLAs and the chess piece. In the particular configuration for 
which Fig.|3] was calculated, this apparent stretching happens as follows. As the chess 
piece is moved further and further away from the CLAs, it first appears stretched in 
the X direction, whereby the x stretch factor changes sign. In fact, this sign change 
occurs at z/e = 5, where the x stretch factor diverges: it becomes +00 as z/e = 5 is 
approached from below, and —00 as it is approached from above. Subsequently the 
chess piece appears very noticeably stretched in the y direction, whereby the y stretch 
factor also changes sign (at Zk = 20). 
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Figure 4. Examples of simulated views through CLAs with different 
combinations of 5x and 5y. The views were calculated with the same parameters 
as those used in Fig. |2] but with rjx = —1/4. 

As the second step of the generahsation we offset the optical axis of the second 
lens with respect to that of the first. (The optical axes of the two lenses remain 
parallel.) The offset d in the two-dimensional case generalizes into offsets dx and dy in 
the X and y direction, respectively. Again, we describe these offsets by dimensionless 
quantities 

4 = ^, 5y = iL, (9) 

Jxl Jyl 

Fig. |4] shows views through CLAs with different combinations of 5^ and 5y. 

As the final generalisation step we rotate the combination of two elliptical lenses. 
We describe this rotation by Euler angles ((/), ^, ?/^), using the most common definition 
of the Euler angles (the x-convention) [30 , in which the general rotation is performed 
in the form of three rotations: first by an angle (j) about the z axis, then by an angle 
9 about the x axis, and finally by an angle ijj about the z axis again. In the following 
we briefly discuss each of these angles, and examples of the effect they have on the 
view through CLAs. 

The angle (j) rotates the major and minor axes of the elliptical lenses away from 
the X and y directions, and therefore provides the generalisation mentioned above. 
Fig. |5] demonstrates the effect of increasing (j) from 0° to 180°. The effect on the 
appearance of the chess piece a finite distance behind the CLAs is a distortion and 
a rotation through 360°. In contrast, the effect on the appearance of the horizon is 
not a full rotation, but a seesawing, as can be seen from the movie in the supporting 
online material (see caption of Fig. [5|. 

The angle 6 describes the inclination of the lenses relative to the planes of the 
corresponding lenslet arrays. It is the 3D generalisation of the angle of the same name 
we discussed in the 2D case (section [sj particularly Fig. [T]). Fig. |6] shows example 
views through CLAs that differ only in their value of 0. Changing shifts the apparent 
position of the chess piece seen through the CLAs in each frame sideways. This effect 
is very similar to that of changing 5x and 5y] specifically, in the case = 0, changing 
6 is approximately the same as changing 5x- This can be seen very clearly in the 2D 
case: if equation ([6| is written for small angles ai, 0^2 and ^, it becomes 

ai-e^3^r]{a2-e), (10) 
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Figure 5. Examples of the effect of a change in the angle (j) on the view 
through CLAs. The frames are taken from a movie showing the view through 
the CLAs as its angle (j) is increased from 0° to 180°. The movie is available in 
the supporting online material (MPEG-4, 320kB). The views were simulated for 
the same parameters as in Fig. 12] but with rjx = —1/4. (The optical-axis offsets 
are zero, i.e. 5x = ^ and 5y = 0, and so are the other angles, i.e. ^ = 0, -0 = 0.) 




e=o° e=o.2° e=o.4° 




e=o.6° e=o.8° e=r 

Figure 6. Examples of the effect of changes in the angle on the view through 
CLAs. The views shown here were simulated with the parameters used in Fig. [s] 
but with (f) = 45° and with different values of as indicated below each frame. 
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Figure 7. View through CLAs with various values of -0. All views were simulated 
for the same choice of parameters as that used in Fig. [6] and with 9 = 0.5°. The 
frames are taken from a movie showing the view through CLAs as its value of 
-0 is increased from 0° to 360°. The movie is contained in the supporting online 
material (MPEG-4, 556kB). 



which can be written as 

Qfi {r]-l)e-5 

Q^2 ^ \ • (11) 

For smah angles, 5 and are therefore not independent but simply act as a combined 
angle that gets added to the scaled angle of incidence, ai/r?, to give the angle of 
refraction, 0^2. 



7| illustrates the effect 
the view depends on 



The angle is required to make the rotation general. Fig. 
of changing this angle. As was already demonstrated in Fig. [3, 
the distance between sheet and object. Fig. |8] demonstrates this again, this time for 
more complex CLAs. 



5. Conclusions 



In this paper we have described generalisations of CLAs [11 that retain their 
important property of acting like a homogeneous light-bending sheet. We have not yet 
investigated in detail the optics of such generalized CLAs, for example their imaging 
properties, but intend to do this in future. 

It is also possible to use the additional parameters locally describing CLAs as 
additional degrees of freedom in general optical-design applications. For example, we 
believe it will be possible to realize optical cloaking [S) [6] from a range of directions 
with structures consisting of CLAs. 

Finally, we intend to build generalized CLAs. 
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z,=16 z,=60 



Figure 8. Effect of changing the distance between the CLAs and the object, 
calculated for more complex CLAs than that for which Fig. [s] was calculated. 
Specifically, the object is a chess piece and its distance behind the CLAs, z^, is 
indicated below each frame. In the final frame (2;^ = 60), the position of the chess 
piece (which appears very small) is marked by a circle. The parameters of the 
CLAs are those of the sheets for which Fig. [T] was calculated, with the additional 
choice if) = 90°. A movie showing the view as is increased from —2 to +38 is 
contained in the supporting online material (MPEG-4, 216kB). The ray-tracing 
image of CLAs shown in Fig. [ijb) was calculated for the same values of rjx^ r]y, 
6x: Sy, (f), 0, and ip, but for a significantly reduced separation between the lenslet 
arrays and correspondingly significantly thicker lenses. 
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